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The amplitude of the magnetic field near the Galactic Centre has been uncertain by two 
orders of magnitude for several decades. On a scale of ~100 pc fields of ~1000 fiG ^ ^ have 
been reported, implying a magnetic energy density more than 10,000 times stronger than 
typical for the Galaxy. Alternatively, the assumption of pressure equilibrium between the 
various phases of the Galactic Centre interstellar medium (including turbulent molecular 
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gas; the contesteci^ "very hot" plasma; and the magnetic field) suggests fields of ~100 fiG 
over ~400 pc size scales^. Finally, assuming equipartition, fields of only ~6 fiG have been 
inferred from radio observations* for 400 pc scales. Here we report a compilation of previous 
data that reveals a down-break in the region's non-thermal radio spectrum (attributable to 
a transition from bremsstrahlung to synchrotron cooling of the in situ cosmic-ray electron 
population). We show that the spectral break requires that the Galactic Centre field be at 
least ~50 /iG on 400 pc scales, lest the synchrotron-emitting electrons produce too much 
7-ray emission given existing constraints'. Other considerations support a field of 100 /xG, 
implying that > 10% of the Galaxy's magnetic energy is contained in only < 0.05% of its 
volume. 

Pinning down the large-scale magnetic field of the Galactic Centre (GC)- which provides our 
closest view of a galactic nucleus - within the greater than two orders of magnitude uncertainty 
implied by existing, rival analyses has long been of interest. A 1000 /iG field would have substan- 
tial ramifications for the region's dynamics including Umiting the diffusion distances of relativistic 
particles (thereby excluding scenarios where the diffuse, ~TeV 7-ray glow of the GCPis ultimately 
explained as due to a single, astrophysical acceleratoi^^. It would also generate a large enough 
magnetic drag to enhance the in-spiral rate of giant molecular gas clouds towards the GC^, limiting 
the lifetime of these to ~100 million years, implying 'star-bursts' on the same timescale. On the 
other hand, the formation of individual stars might be inhibited (or the stellar initial mass function 
biased) by a strong magnetic field's tendency to support gas against gravitational collapse. 
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Radio observations at 74 MHz and 330 MH2P reveal a diffuse (but distinct) region of non- 
thermal radio emission covering the GC (out to ~ ±3° or ~ ±420 pc from the GC along the 
Galactic plane). Invoking the "equipartition" condition (minimizing the total energy in magnetic 
field and relativistic electrons), a field of only 6 is inferre(i^, typical for the Galaxy-at-large, 
climbing to 11 fiG over the inner ~ ±0.°8 ( for extremal parameter values the field might reach 
100 /iG). 

To probe this radio structure at higher frequencies, we have assembled total-intensity, single- 
dish flux density data at 1.4, 2.4, 2.7, and 10 GHz?^^. These data are polluted by line-of-sight 
synchrotron emission in the Galactic plane both behind and in front of the GC and the flux density 
contributed by discrete sources; we used a combination of low-pass (spatial wavenumber) filtering 
to remove the latter and aU-sky Galactic synchrotron background observations to remove the for- 
mer (see Supplementary Information). After this processing a distinct, non-thermal, radio structure 



is revealed in all radio maps up to 10 GHz (see Figure 1 1. However, the structure's 74 MHz-10 



GHz spectrum is not described by a pure power-law: attempting to fit such to the cleaned data we 
find a minimum of 4.9 per degree of freedom (dof = 4), excluded at a confidence level of 3Aa 
( [Figure 3[ ). 



In fact, fitting separate power laws to, respectively, the background-subtracted lower three 
and upper three radio data, we find that these extrapolations intersect at ~ 1.7 GHz with a spectral 
break of ~0.6. This is close to the canonical break of 1/2 produced by a steady-state synchrotron- 
radiating electron population that transitions (with increasing energy) from bremsstrahlung-cooled 
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to synchrotron-cooled. Since the synchrotron coohng rate is a function of magnetic field, B, and 
bremsstrahlung a function of uh, the break frequency is a function of both parameters. Observation 
of a spectral break then determines acceptable pairs of B and hh for the synchrotron-emitting 
electrons' environment. 

We have modelled the cooled electron distribution and resulting synchrotron emission as a 
function of B, uh, and spectral index at injection of the electron population. We have accounted 
for losses due to inverse-Compton emission following collisions with ambient lighP^, ionization, 
bremsstrahlung and synchrotron emission. 



The results of our fitting procedure are displayed in Figure 2 and Figure 3 A very important 
constraint is proffered by 7-ray data covering the non-thermal emission region: bremsstrahlung and 
inverse Compton emission will inescapably be generated by the electrons responsible for the ob- 
served radio emission (the energies of electrons synchrotron-radiating at ~GHz and bremsstrahlung- 
radiating at > 100 MeV are very similar). This must not surpass the 300 MeV 7-ray flux from the 
region measured by EGREl'^ (see Figure 2 This consideration rules out field amplitudes < 50 
/iG. 

In the context of diffuse, ~GeV emission around the GC, we eagerly await results from the 
Fermi Gamma-ray Space TelescopeP' whose sensitivity (more than an order of magnitude better 
than egret's) promises, at worst, an increased lower-limit to the large-scale magnetic field or, 
more optimistically, a measurement of this field (in concert with radio observations). For now we 
emphasis that it is a novel analysis - not new data - that has allowed the new constraint on the GC 



magnetic field. 



Figure 4 shows tlie energy densities of various GC ISM phases. For acceptable field am- 
plitudes, the cosmic ray electron population is considerably sub-equipartition with respect to the 
other GC ISM phases (even after accounting for filling factor effects), in particular, the magnetic 
field, explaining why the magnetic field estimate arrived at assuming equipartition^ is too low. We 
have also calculated the maximum possible energy density in the GC cosmic ray proton population 
given the EGRET 7-ray constraints. Note that at ~ 100 yuG the magnetic field reaches equiparti- 
tion (at ~300 eV cm"^) with the putative 'very hot' (~ 8 keV) phase of the X-ray-emitting plasma 
supposedly detected throughout the central few degrees along the Galactic planeP^. Moreover, the 
energy density of the gas turbulence kinetic energy for the derived uh is within a factor of a few 
of equipartition with these other phases up to magnetic field amplitudes of ~ 100 fxG. 

This situation - implying near pressure equilibrium between a number of GC interstellar 
medium phases (including the plasma) - mirrors that of the Galactic disk, albeit at much higher 
pressure. Such considerations led to the predictioiP that the real GC magnetic field amplitude lies 
close to ~ 100 /iG. Very recently, however, observations have shown"^ that the X-ray emission 
from a region around / = 0.08, b =1.42 (taken to be typical of the so-called X-ray Ridge) is due to 
unresolved point sources, implying the very hot X-ray plasma is illusory. This casts the pressure 
equilibrium argument^ into some doubt. Note, however, that our magnetic field lower limit holds 
irrespective of whether the 8 keV plasma is real or not. 

There are a number of intriguing parallels above to the situation apparently pertaining within 
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the inner regions of starburst galaxies. These independently support the idea that the field is ac- 
tually close to 100 fiG. In starburst environments, it is contended^, equipartition magnetic field 
values obtained from radio observations significantly underestimate the real field. Fields are actu- 
ally sufficiently high to be in hydrostatic equilibrium with the self-gravity of the gaseous disk of the 
starburst. Such strong fields (together with high gas densities) guarantee that the relativistic parti- 
cle population dumps all its energy before being transported out of the system. This calorimetric 
hmip2 may explaiiP why even very luminous star-bursting galaxies fall on the far-infrared-radio- 
continuum correlation^^. Circumstantial evidence also suggests that radio emission from starbursts 
is dominated by secondary electrons created in collisions between cosmic ray ions and gad^ (rather 
than directly-accelerated electrons). 

It is noteworthy, then, that the 60 micron and 1.4 GHz emission from the radio emission 
region place it within scatter of the far-infrared-radio-continuum correlationP^. Furthermore, in 
the GC a ~100 fxG field is precisely in the range required to establish hydrostatic equilibrium 
(given the total and gaseous surface densities). A situation wherein a magnetic field provides 
significant pressure support against gravity may lead to the development of the Parker instability^ 
and exactly this is suggested by mm- wave observation^ of molecular filaments of several hundred 
parsec length within ~ 1 kpc of the GC. These observations independently suggest a ~ 100 /iG 
field. Finally, the diffuse, ~TeV 7-ray glow from the vicinity of the GCPis most likely explained by 
cosmic ray impacts with gas. Unavoidably, such collisions would also produce copious secondary 
electrons which could then contribute significantly to the region's synchrotron radio emission (for 



fields > 300/iG, 100% of the radio emission from the rectangular region shown in Figure 1 could 



be attributed to secondary electrons). Taken altogether, these facts paint the Galactic centre as akin 
to a weak starburst with a magnetic field of ~ 100 yuG. 
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Figure 1 | Total intensity image of the region at 10 GHz. Radio map convolved to a resolution 
of r.2 X 1°.2 with contours at 10, 20, 40, 80, 160 and 240 Jy/beam. (Native resolution and 
convolved images at z/ > 1.4 GHz in the Supplementary Information.) There is a striking constancy 
in the appearance of the radio structure from 74 MHz to at least 10 GHz (the large ellipse traces 
the diffuse, non-thermal radio emission region first identified at 74 and 330 MHz^). The small 
rectangle delineates the region from which the HESS collaboration determines a diffuse ~ TeV 
7-ray intensit)^. 
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Figure 2 



Figure 2 | Spectrum of the region: data and models, (a) Radio. Points:-circles: flux den- 
sity of the radio structure after removal of flux from sources with sizes < 1°.2. The vertical 
line divides the 74 MHz datum (that receives no contribution from Galactic plane synchrotron 
background/foreground) from all the other data (that do) - hence the discontinuity in the model 
curves at 100 MHz; squares: discrete source flux densities (measured directly by the VLA at 
330 MHzP] and otherwise obtained through the Fourier technique described in the Supplemen- 
tary Information); arrows: discrete source lower limits obtained from ATCA. The error bars show 
68% confidence intervals. Curves: the best-fit flux density due to synchrotron emission from a 
cooled electron distribution plus Galactic plane synchrotron for the case of (solid) 100 fiG and 
(dotted) 30 (j,G fields (the break in the radio synchrotron curve mirrors a corresponding break in 
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the distribution of radiating electrons; we assume the electrons are injected with a power law in 
momenturnp3); dashed: the best-fit null hypothesis case of a pure power-law signal plus Galactic 
plane synchrotron, (b) Broadband. Models for magnetic field of 100 /iG (solid) and 30 /xG (dot- 
ted). The injected electron spectrum is assumed to follow a power law up to ~ 100 TeV. (Note 
that the cut-off energy is only weakly constrained. Our conclusions are not sensitive to the exact 
cut-off energy, however, given that radio observations guarantee the existence of ~GeV electrons 
whose bremsstrahlung emission we can compare against the EGRET intensity datum.) The lower 
energy curves are synchrotron emission. The higher are bremsstrahlung plus inverse Compton 
emission from the same electrons responsible for the synchrotron (a figure showing the individual 
bremsstrahlung and inverse Compton contributions is shown in the Supplementary Information). 
The upper limits are due to observations by EGREl^i at 300 MeV and HES^ at TeV. 
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Figure 3 | Plot of the per degree of freedom as a function of magnetic field amplitude. 

Curves: (as labelled) model for a cooled primary electron model (with 3 degrees of freedom) and 
the null hypothesis of a pure power law electron distribution (with 4 dof). The solid curves are 
constrained by the requirement that the 7-ray emission from the synchrotron-radiating electron 
population be less than the upper limit obtained from EGRET data' (dotted curves not so con- 
strained). The horizontal dashed lines mark the l,2,3,4cr confidence limits for a model with 3 dof 
and do not apply to the null hypothesis (which only achieves a best fit acceptable at the ~ 3.4cr 
level). 
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Figure 4 | Energy density in various phases of the GC interstellar medium as a function 
of magnetic field. Bands show 68% confidence limits (except for p). Labels denote:- 'mag': 
magnetic field; 'plasma': the disputed^ X-ray emitting 8 keV plasma (with a density*^ 0.03-0.06 
cm"^); 'turb': turbulent motions of the local gas (assuming it is at the best-fit rin and has a veloc- 
ity dispersion in the range 15-30 km/s typical for GC molecular cloud^- note that, as the uh 
required to fit the radio spectrum increases with B, the turbulent energy density is an increasing 
function of B)\ 'e': cosmic ray electrons; and 'p': a conservative upper bound on the cosmic ray 
proton energy density inferred from the 300 MeV 7-ray upper limit (a putative proton population, 
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colliding with ambient gas at the best-fit uh, would produce 7-rays mostly via neutral meson de- 
cay). The electron and proton energy densities take into account the filling factor of gas at the 
given value of nn^i (we take the filhng factor of gas within la around the best-fit uh value); 
this is a correction upwards to the values of and Up by an amount 100-1000. Without this 
correction, one would have f/g — f^B at ~ 10 /iG in agreement with the estimate derived assuming 
equipartion^. A factor of 2 uncertainty^^ in the hydrogen mass of the region is also accounted for 
in determining the electron and proton energy densities. 
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